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Abstract 

It is known that, in the chiral Umit, spatially inhomogeneous chiral condensate 
occurs in the Nambu-Jona-Lasinio (NJL) model at finite density within a mean- 
field approximation. We study here how an introduction of current quark mass 
affects the ground state with the spatially inhomogeneous chiral condensate. Nu- 
merical calculations show that, even if the current quark mass is introduced, the 
spatially inhomogeneous chiral condensate can take place. In order to obtain 
the ground state, the thermodynamic potential is calculated with a mean-field 
approximation. The influence of finite current mass on the thermodynamic po- 
tential consists of following two parts. One is a part coming from the field energy 
of the condensate, which favors inhomogeneous chiral condensate. The other is 
a part coming from the Dirac sea and the Fermi sea, which favors homogeneous 
chiral condensate. We also find that when the spatially inhomogeneous chiral 
condensate occurs, the baryon number density becomes spatially inhomogeneous. 
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1 Introduction 



Quantum chromodynamics (QCD) at finite density is one of interesting topics in tliese 
days. Tlie study of tliis field will help us to understand the physics of neutron stars, 
compact stars, and heavy ion collisions. In the vacuum state of QCD (quark chemical 
potential = 0) at low temperature, the chiral symmetry is broken spontaneously 
and the confinement occurs. On the other hand, at extremely high densities {fi has 
very large values) and low temperature, color superconductivity will be realized in QCD 
[H [2] . At such high density region, the perturbative calculation with the gauge coupling 
g is possible due to asymptotic freedom. 

Now, at the moderate density region (moderate value of fi) where the coupling g 
is not small, how a ground state of QCD becomes? The method of perturbation of the 
coupling g cannot be used, and also it is difficult to apply the lattice QCD simulations to 
a system with finite density. Therefore, people have used effective theories of QCD such 
as the Nambu-Jona-Lasinio (NJL) model [21 S] to study the physics at the moderate 
density region[5]. One of the interesting topics of these researches is that spatially 
inhomogeneous chiral condensate occurs in the ground state at the moderate density 
pin [151 [20!. 

Before discussing this topic in detail, let us look back the idea that fermionic con- 
densate becomes spatially inhomogeneous at finite density. The possibility of spatially 
inhomogeneous chiral condensate in QCD at finite density was first discussed in RefjSj. 
The authors in Ref[H] have shown that in the limit, N^. (number of colors) —>■ oo, the spa- 
tially inhomogeneous chiral condensate occurs in the ground state of QCD at extremely 
high density region such that the relation g'^Nc <^ 1 holds. That state is the standing 
wave ground state having the wave number 2/i, in which particle and hole with the 
same Fermi momentum p{\p\ = /u) condense (Overhauser effect [H]). When one takes 
the number of colors A^^^ of QCD to the realistic number three, however, it is shown 
by the later researches that the BCS effect (particle-particle condensation) is superior 
to the Overhauser effect (particle-hole condensation) and the color superconductor is 
realized at extremely high density [TOl [Til [I2j. In some two-dimensional models, it is 
argued that spatially inhomogeneous chiral condensate also occurs at finite density. In 
the chiral Gross-Neveu model with the limit ^ oo at finite density, it is shown that 
both scalar {ibip) cind pseudoscalar {ipi 75 ip) become spatially inhomogeneous in the 
ground state p| [T3l [11]. 

Now, let us return to the topic mentioned above. In Ref [HI [7], the NJL model with 
Ac = 3 and Nf (number of fiavors) = 2 is considered at moderate density. Assuming 
the following mean-field, 

(ipip) = A ■ cos{q ■ r), {ipi-f^Tstp) = A ■ sin {q ■ r), (1) 

^ In the chiral Gross-Neveu model, finite bare quark mass case has also been studied [13]. When 
the bare quark mass is taken to be finite, the chiral angle depends on the value of the bare mass and 
the baryon density becomes spatially inhomogeneous. 
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where g is a wave number vector, the authors in Ref P, [7] obtain the ground state 
of that model by finding the minimum value of the thermodynamic potential in the 
mean-field approximation. In the chiral limit, they find numerically that the ground 
state has non-zero value of q at low temperature and a high density region, namely 
spatially inhomogeneous dual chiral condensate, Eq.([T]), occurs at that density region. 
Furthermore, the quark number density (ip^ip) becomes spatially homogeneous in the 
chiral limit. 

In Ref [20] , the T — phase diagram of the N JL-type model is studied around the 
chiral critical point in the chiral limit, which includes spatially inhomogeneous chiral 
condensate. There, the condensate is restricted to only scalar form instead of allowing 
both scalar and pseudoscalar form like Eq.([T]). With more general ansatz for the scalar 
condensate, two phase transitions are found |20]. One is from the homogeneous mas- 
sive phase to inhomogeneous phase, and the other is from the inhomogeneous phase 
to chirally symmetric phase, both of which are the second order. In the real world, 
however, the current quark mass is not zero and the pion (NG boson) has a mass of 
about 135MeV. Recently, the same issue is studied for the finite current quark mass 
|15j . In general, the thermodynamic potential is a function of the complex order 
parameter M(r), 

M{r) = m-2Gs{{tp4j) + i{4Jii5r3^)}, (2) 

and the stationary constraint, SQ/SM{r)* = 0, is held in the ground state. In Ref |15j . 
the stationary constraint is solved when the order parameter M(r) is restricted to real 
(i.e., pseudoscalar condensate = 0). The numerical calculation using that real solution 
shows that the qualitative feature of the phase diagram remains unchanged. Now, how 
about the case where the order parameter M(r) takes a complex value (i.e., both scalar 
and pseudoscalar condensate)? Unfortunately, it is very difficult to solve the stationary 
constraint with finite current mass when M(r) is a complex value. Although no solution 
is known, we need to study the case where the order parameter M(r) is complex by, 
for example, assuming ansatz for M(r). 

In this paper, we study the massive NJL model in the case where the order pa- 
rameter is complex, that is, both scalar and pseudoscalar condensate exist. Concretely, 
we study how an introduction of finite current mass affects the ground state with the 
spatially inhomogeneous dual chiral condensate, Eq.([T]), in Ref [HI [7|. Does the dual 
chiral condensate remain spatially inhomogeneous? And, does the quark number den- 
sity remain spatially homogeneous? To find out these problems at moderate baryon 
density, we use the NJL model with finite current quark mass [15] (not the chiral limit 
case as in Ref [6l[7|). Here, let us compare the ansatz, Eq.([T]), and that in Ref [20]. The 
ansatz, Eq.([T]), allows the pseudoscalar condensate to occur in addition to the scalar 
condensate, while the ansatz in Ref [2i)| restricts the condensate to the scalar form. In 
this respect, we can say that Eq.([T]) is more general than the ansatz in Ref [20]. On the 
other hand, concerning Fourier analysis of the assumed condensate, Eq.([T]) is restricted 
to monochromatic wave, while the ansatz in Ref [20] contains various higher harmonics. 
From the viewpoint of this, we can say that the ansatz in Ref [20j is more general than 
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Eq.([T]). Now, we return to the discussion of the NJL model. In the chiral hmit, the 
thermodynamic potential can be calculated P, [7] in a mean-field approximation under 
the assumption Eq.([T]). On the other hand, when the current mass is taken to be finite, 
it is difficult to calculate the thermodynamic potential because the fermion propagator 
depends on the space coordinates r explicitly. In order to avoid this difficulty, we will 
expand in powers of the current mass m. The thermodynamic potential uj will be cal- 
culated up to first order 0{m) and the ground state is obtained by finding a minimum 
value of u. 

This paper is organized as follows. The NJL model with Nf. = 3 and Nj = 2 is 
introduced in section 2. We assume that mean fields (ipip) and {ipi 75 ip), Eq.([T]), exist 
at finite density, and the thermodynamic potential is calculated analytically in a mean- 
field approximation. In section 3, the numerical calculations of the thermodynamic 
potential at zero temperature are carried out and we obtain the ground state which 
minimize the thermodynamic potential. For a given quark chemical potential fi, these 
numerical calculations enable us to find whether the spatially inhomogeneous chiral 
condensate is realized {q 7^ 0) or not (q = 0). In section 4, we show that the quark 
number density becomes spatially inhomogeneous when the spatially inhomogeneous 
chiral condensate is realized in the case of finite current quark mass. Section 5 is 
devoted to conclusions. 



2 The Nambu-Jona-Lasinio model and 
standing wave ansatz 

The Lagrangian of the NJL model is 

C = - m)^ + G + (T/^i 75 Ti)f\ , (3) 

where the number of colors is iV^ = 3 and the number of fiavors is N ^ = 2. Note that 
the current quark mass is taken to be finite and we set = rrid = m > 0. Now, we 
suppose the following mean fields exist at finite density [H [7] , 

{^{J^P) = C cos (q-r), {^/ji ^) = C sin {q ■ r) , (C < 0), (4) 

where q is the wave number vector and C is a constant. The other components are 
assumed to vanish, {ipi 75 ri ip) = {ipi 75 T2 ip) = 0. When q vanishes, these have the 
usual forms (ipip) = C and {tpi 75 ip) = 0. These mean fields are put on the chiral 
circle, (ipip)'^ + (V'«75 T3 ip)"^ = C^. 

Here we shall make a comment on the ansatz Eq.(jl]). If the chiral limit, m — 0, 
is considered, the bottom of the effective potential is the chiral circle, and it would 
be natural for (ipip) and {ipi 75 ip) to be put on that chiral circle. However, when 
the current mass m is finite, the effective potential slightly tilts in the direction of 
't/jip and the shape of the effective potential becomes very complicated. In such case. 
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there might exists another ansatz which has more appropriate function of r than the 
trigonometrical function used in Eq.(jl]). This point will be discussed in subsection 3.1 
and section 5. In the present paper, we use the ansatz Eq.(jl]) for simplicity. 

One can choose the direction of the wave number vector q as q = (0, 0, q),{q > 0) 
without loss of generality. The mean-field approximated Lagrangian becomes 

'Cmf = ^[2^-m + /i7o-Mexp(z75r3q-r)]^- — , (5) 
where the parameter M has been defined as 

M = -2GC = 2G^{^^y + {^t^5r3ijy, (6) 

and /i is the quark chemical potential. In the path integral representation, the bilinear 
form of fermion in Cmf can be written as 



V ijjVi/j exp i J d xip [i<^ — m + n'fo — M exp (275 r^q ■ r)] ijj 
J V ip'Vip' exp i J d^x i)' i-f^'id^ + ^75 r^q^) - M + /i 70 - m exp (-^75 r^q ■ r) 



X exp i J d^x iIj' 27^(9^ + ^75 r^q^) - Mt + fi-fo- m {exp (-^75 r^q ■ r) - 1} 



(7) 



where '0' = exp |+| 75 Tsq-rj^/' and M + m = Mt > m. The propagator 5* of the field 
ip' is 

g _ , }i /gN 

il^id^ + f 75 T-sg^) - Mt + /i 7o - m {exp (-^75 r^q ■ r) - 1} ' 

If the wave number vector q vanishes, S represents a propagator of free fermion with 
a mass M^. 

When the wave number vector q 7^ and the current quark mass m 7^ 0, the 
denominator of the propagator S depends on the space coordinates r explicitly, so that 
it becomes difficult to calculate the expression involving S. To avoid this difficulty, we 
separate 5*"^ into two, 

= So' - Vm, (9) 



where 



Using the identity 



Vm = m{exp{-i-f5T3q ■ r) - 1} . (10) 
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for any noncommutable operators A and B, we obtain 



S = So + SoV„,S. (12) 

From this equation, we have 

S = Sq + SoVmSo + SoVmSoVmSo + • • • . (13) 

Here, it should be noted that in the chiral hmit, m ^ 0, we have, — > 0. Since 
the current quark mass m is the smallest energy scale in the system, we treat the 
Vm ~ 0(m) as perturbative part and neglect higher order O(V^). We shall use the 
approximate equation (fT3!) to calculate the expression involving S, as will be done 
in section 4. For the expression involving S~^, the separation Eq.([9]) will be used as 
follows. 

The thermodynamic potential u is useful so as to find the ground state of the 
system with finite density fi > and finite temperature T > 0. In the imaginary time 
formulation [TB], r = it, the mean-field approximated thermodynamic potential uj has 
the form, 

T 

oj = — 



— \nN' J V ipVip exp J J d^x Cuf 

— In / ViP'ViP'exp / dr / d^xip'S-^ip'. (14) 

AG V J Jo J 



Here V = is the volume, f3 is 1/T, and the irrelevant constant is omitted. Since 
the path integral calculation in Eq. (fT^ is difficult, we regard the term Vm in Eq.([9]) 
as perturbative part and obtain u up to the order 0(Kn). We separate the action 
Jq dr J d^xijj'S'^ip' = A into two parts referring to Eq.([n]), 



A = Ao + Am. (15) 



Here, Aq and Am are 



Ao = J^^ dr J d^x^'S^^ip', 

Am = j^dr J d'x^'{-Vm)^', (16) 

and note that in the chiral limit, m —>■ 0, one has. Am 0. Expanding exp{Am), we 
have 



w 



In y D ijj'Vip' exp A 

--ln/D^'P^'exp(Ao)E77^!n 

1 /■-n/'-n/' , , JV^'V^'Amexp{Ao) , 



M2 


T 


4G 


~ V 




T 


4G 


" V 


M2 


T 


4G 


" V 



(17) 
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The leading term in Eg. (11 71) is 

Inj V^'V^'exp{Ao) 



--NJnDet 

V It 



= -TNr 



= UJB+ UJf, 



where 



1 



{p^ = iujj + fx = i{2j + l)7cT + fi). (19) 



° - Mt- n ' 

The propagator Sq has four energy poles, po = e„(ri = 1, 2, 3, 4), with 



ei = e_ = 

€2 = e+ = 

es = -ei, 
£4 = -62 • 



+ >o, 



IPl+{\fMf+?. + l 



> 0, 



The ub and t^p have the form P, [7j , 



ujB = -NcNfjj^ (ei + e2) + rin|l + exp 



(27r)3 



T In 1 + exp 



+Tln 1^1 + exp 



T / 

(g2 + /i) 

T 



T 



+T In 1 + exp 



(£2 - /i) 
T 



(20) 



(21) 



(22) 



The ujb represents the contribution from the Dirac sea and the contribution from 
the Fermi sea. The next leading order of 0{Vm ) is [I6] 



= -^iVe(-l)m 
= 5u;d + Sujf, 



i MLq/2) Ayy^f d'P . . 



(23) 
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where [I7j 
6ujb = 



'sm{Lq/2) 
. {Lq/2) 



1 1 m N^Nf j 



(2^ 



X 



E 

n=3,4 



Mt{JMf+pl 



■sin(Lg/2) 
. {Lq/2) 



1 + exp 



(en - 



(24) 



6uj-p = 



(Lg/2) j ^ ^7 (27r)3 



X 



E 



n=l,2 

■sin(Lg/2) 



d 



I (^?/2) 

Both 5ujy) and ^cjf vanish in the chiral hmit, 

hm 5uj-D = hm 6ujf = 0, 

m— >0 m^O 

as well as in the limit, g — 0, as it should be. 

Eventually, thermodynamic potential uj up to 0{Vm) becomes 



UJ = UJD + UJp + 



AG 



m/1 



-1 



(25) 



(26) 



(27) 



where 5ujd and 5uj-p come from the term, Vm- In the zero temperature limit, T ^ 0, 
each term becomes 



^D = -iVciV//^(6i + e2), 

a;F = -iVciV/ / 7^ [ (/i - ei) 0(/x - e,) + (/i - 62) 0(/i - e^)] 



sin(Lg/2) 1 9 
— ^ — 1 > m - — 



(Lg/2) 
■sin(Lg/2) 
. (^g/2) 



9Mi 
d 



(28) 
(29) 
(30) 
(31) 



Note that cuf is always negative, u-p < 0. But 5ujy can have either positive or negative 
value. 
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3 Numerical calculation with zero temperature 



When the values of chemical potential /i and temperature T are given, the ground 
state of our system can be described by two parameters (Mf,g) which minimize the 
thermodynamic potential uj. We will obtain the ground state in the zero temperature 
case T = by numerical calculations. Because the NJL model is a cut-off theory, 
a regularization method should be specified to define the theory, and we utilize the 
proper-time regularization method here. With the proper-time regularization, uj-q 
has the following form P [7] , 



°° dk, f°° dr 



i dkz f 



1/A2 



exp <; - ( VA:2 + ^ 



q ^ ^ 



r 







M 


II 







(32) 



where A is a cut-off parameter. The regularized Sud can be obtained from Eq. (l32l) by 
use of Eq.(|2il). 

We now set the three input parameters {G, m, A) so that the observed values /,r = 
92.4 MeV and m^r = 135 MeV are reproduced in the vacuum state /i = 0. We first give 
arbitrary values of A, and then fix G and m so as to suit /jr and m,r to their observed 
values. In Appendix A, the way of determining {G, m, A) is given in detail. There still 
remains one degree of freedom, because we use only two observable quantities /^r and 

to determine the three input parameters. Therefore, we require further that the 
chiral phase transition should be first order in spatially homogeneous case {q = 0) at 
finite density. The parameters {G, m, A) obtained in Appendix A are restricted by this 
requirement. The spatially homogeneous case [q = 0) will be considered in the next 
subsection. 



3.1 Spatially homogeneous case at finite density 

Before studying spatially inhomogeneous case {q 0), we consider in this subsection 
spatially homogeneous case [q = 0), which is easier to deal with. As is well known 
|19j . chiral symmetry is restored as the quark chemical potential fi increases (density 
increases). Here we restrict the values {G,m,A) obtained in Appendix A by requiring 
the chiral phase transition to be first order when g = 0. For a given value of fi 
(and T = 0), the ground state is described by the parameter Mt which minimize the 
thermodynamic potential uj. By the numerical calculations, we find that the chiral 
phase transition becomes first order when the input parameters {G, m, A) obtained in 
Appendix A satisfy the following relation, 

GA^ > 5.6. (33) 

Henceforth, we use the input parameters' values (G, m. A) obtained in Appendix A 
which also satisfy the condition Eq. (l33|) . 
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In analyzing the behavior of Mj, we find the following by numerical calculations. As 
the chemical potential fi increases, the parameter Mj decreases discontinuously and the 
first order phase transition of chiral symmetry occurs. If fi increases further (/i < A), 
one might expect that the parameter Mt will decrease more and approach m, Mt > m. 
However, our observation is that, after chiral phase transition, Mt decreases more and 
eventually it becomes Mt < m when exceeds about 0.7A. The relation Mt < m 
implies that the dynamically obtained quark mass Mt is less than the current quark 
mass if /i ~ 0.7A, and this phenomenon is unnatural. Hence we restrict ourselves to 
the region, /i ~ 0.7A and do not consider the value of /i larger than 0.7A. The fact 
that the relation Mt < m holds in the region /i ~ 0.7A will be caused by our choice of 
the regularization method, that is the proper-time regularization. This point will be 
discussed in Appendix B. 

Before closing this subsection, we would like to study the effects of current mass 
on the spatially homogeneous chiral condensate. This search would be helpful to in- 
vestigate the effects of current mass on the spatially inhomogeneous chiral condensate 
which will be discussed in section 3.3. In the ansatz, Eq.(jl]), when the chiral conden- 
sates are spatially homogeneous (g = 0), the condensates are chosen to be {il)ip) = C 
and {ipi'j5T3ip) = 0. In the spatially homogeneous case, it is known that this choice 
realizes the lowest value of the energy if the current quark mass is finite. Here, we make 
sure of this fact by the thermodynamic potential uj expanded in Vm- To this end, we 
assume in this subsection the following spatially homogeneous chiral condensate, 

(^^) = C cos 0, {ipi 75 T^'^) = C sin 0, {C <0), (34) 

where the chiral angle is a constant ( < < 27r ). The thermodynamic potential 
having these condensates can be calculated as done in Section 2, 

+ + ^1 + {-i^Mt + 5u^ + 5ij^ + — \+ 0(0, (35) 

where Vm = m {exp (—^75 T30) — 1}. Each term has the following form, 

= -N^Nf I ^ 2 Vp2 + M2, (36) 

u;^ = -N,Nf J ^2 (/.-v/p^ + M,^) v/p^ + M^), (37) 

5cud = (cos 0-1) m{-l)N,Nf f 2 ■ ^* 



(27r)3 ^p2 + 



d 

(cos0 — 1) m 7— — tUD, (3J 



dMi 



<5u;f = (cos0-l)mA^eAr; / 7^2 _ g ("^ _ y^p2 + 

d 

= (cos — 1) m cup- (39) 



dMt 
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Thermodynamic potential uj is regarded as a function of two variables, {Mt,(j)), and 
the ground state is described by (Mj, 0) which minimize u. Since the u depends 
on the chiral angle through cos0, one can constraint the chiral angle, < < vr. 
The part, {ujd + up + /AG}, in Eq. fl35l) does not depend on (p. The terms which 
represent the current mass effects are {—{m/2G)Mt + ^cud + Sup} (the term, w? /AG, 
can be regarded as a constant because it does not depend on Mt and 0). The term, 
— {m/2G)Mt, depends on only Mt, and this term lowers u. Next, the term, 6ud, takes 
positive value, Sujd > 0. This 6ujy) is a monotone increasing function of the chiral angle 
in the range of < < vr, therefore the term,(5to'D, tends to move the chiral angle 
closer to the value zero. Thirdly, the term, 6ujf, vanishes if fi < Mt, whereas it takes 
negative value 6ujp < if /i > Mt. This Sup is a monotone decreasing function of in 
the range of < < vr. Therefore the term, Sup, tends to move closer to the value 
TT, and this tendency rises when the chemical potential fi becomes larger. We can say 
that the term Sub and the term 6ujf compete. Then, how about their sum, {6ujf + 6ujy)) 
? The regularized 6ujb is obtained from Eq. fl32l) by use of Eq. fl38l) . and the sum is 



6ljf + Sujb = (cos - 1) 2 m N^NfMt (-1) 

1 „-M?r f 



^7r2 Ji/A2 r2 J (27r)3 1^2 + 



(.40) 



The equation enclosed in the brackets of the right-handed side is the same with that 
in the brackets in Eq. flB.ip . This equation takes positive value when < fi < fic, and 
becomes zero if fi = fic, as is discussed in Appendix B. The equation is a monotone 
decreasing function of /i. Here, fic is defined in Eq.f lB.2p and takes the value, fic ~ 0.7A 
if Mt/A ^ 1. The sum, {Sup + Sud), is then a monotone increasing function of in 
the range of < /i < /ic. We therefore conclude that, in the range of < fi < fic, 
{6ujp + Sub) tends to move the chiral angle closer to the value zero, and this tendency 
becomes weaker when /i grows. From the above, we have shown that, in the range 
of < fj, < fj,c, the condensate, {ipi') = C, {ipi^^T^il:) = 0, realizes the lowest value 
of the thermodynamic potential u if the chiral condensate is restricted to be spatially 
homogeneous. 

The foregoing computation helps us discuss qualitatively the appropriateness of the 
ansatz Eq.(jl]) which we put for the spatially inhomogeneous chiral condensate. When 
H takes the value zero, the above computation concerning the spatially homogeneous 
condensate shows that the u does depend on the chiral angle and one cannot ignore 
that dependence. Then Eq.(jl]) would not be appropriate if /i has small value because 
(ipip) and (ipi'j^T^ip) are put on the chiral circle in Eq.(jl]). On the other hand, the 
dependence of u; on becomes weaker when the value of n becomes larger. Therefore 
we expect that the ansatz putting (tptp) and {ipi'^^T^ip) on the chiral circle such as Eq.(jl]) 
is a moderate assumption when /i takes large value if we neglect the contribution of 
order 0{V^). 
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Figure 1: The mass Mf (solid) and the wave number q (long-dashed) as a function of the 
chemical potential fi, in units of A. For comparison, Mt(q — 0) without inhomogeneous 
chiral condensate is also represented by dotted line. 

3.2 Numerical results 

As is discussed in the previous subsection, we restrict ourselves to the following region, 

GA^ > 5.6, II ~ 0.7A, (41) 

the first condition comes from the requirement that the chiral phase transition is first 
order when q = 0, and the second comes from the requirement that the dynamically 
obtained mass Mt is not less than the current quark mass m when q — 0. The parameter 
L is taken to be a enough large value, L — 2x 10^/A. For a given value of the ground 
state is obtained by determining two values {Mt, q) which minimize the thermodynamic 
potential u. If g = is the solution which minimizes u, the chiral condensate is spatially 
homogeneous. On the other hand, if g 7^ is the solution which minimizes cu, the chiral 
condensate is spatially inhomogeneous. According to numerical calculations, we always 
have the solution q — ior arbitrary values of // when 5.6 < GA^ < 6.94. Namely, in 
the region 5.6 < GA^ < 6.94, the ground state is not spatially inhomogeneous chiral 
condensate. 

If GA^ is larger than 6.94, GA^ > 6.94, however, there exist the solutions q ^ 
for large values of fi. Here we show a concrete example of numerical results with the 
input parameters, GA^ = 7.513, A = 635MeV, and m = 15.3MeV, {Mt (/x = 0) = 
380MeV) in Figure. 1. In the low density region, < ^ < 0.5885A, we have g = and 
the chiral condensate is spatially homogeneous as usual. In the high density region, 
0.5885A < < 0.7A, we get q ^ (long-dashed line) and spatially inhomogeneous dual 
chiral condensate is realized. Let us see the behavior of Mt represented by solid line. 
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Since Mt decreases discontinuously at /i = 0.5885A, the chiral phase transition is first 
order. In this high density region, fi > 0.5885A, one finds by numerical calculations 
that the relation, /i < Mt + q/2, holds. Therefore, the €2 state is occupied with no 
quark because of €2 > Mt + q/2, and all quarks in the Fermi sea are filled in the ei 
states. For comparison, we also show a local minimum solution Mt{q = 0) with q = 0, 
which is represented by dotted line in Figure. 1. In the region, fi > 0.5885A, the points 
(Mj(solid line), q 7^ O(long — dashed line)) in the Mt—q plane correspond to the absolute 
minimum energy state (ground state) of u, and the points (Mt (dotted line), g = 0) 
correspond to the local minimum energy state. The spatially inhomogeneous (g 7^ 0) 
dual chiral condensate occurs at densities higher than about the usual (g = 0) chiral 
phase transition point, /i = 0.59A. The characteristic feature in Figure. 1 is as follows. 
As the chemical potential n increases over the value of 0.5885A, Mt decreases a little. 
Even at high density of = 0.7A, the value of Mt remains rather large, Mt ~ 0.17A. 
This feature is different from the result obtained in the chiral limit case O [?]. If GA^ 
takes the values, 7.04 < GA^ < 8.1, the behavior of q (Mt) in the q — plane (in the 
Mt — ij, plane) is similar to that in Figure. 1. On the other hand, if GA^ takes the values, 
6.94 < GA"^ < 7.04, the wave number q remains zero at the value of chemical potential 
fichi at which Mt decreases discontinuously. The wave number q becomes non-zero at 
the value of the chemical potential whose magnitude is over /Xchi- 



3.3 Effects of current mass 

In this subsection, the effect of current mass on the thermodynamic potential u, 
Eq.(!27I). with T = is cleared. We divide u into two parts, 

a;D + a;F + ^1 + \-l^Mt + (5^D + 5^f) + 4^] + 0{V^). (42) 

The thermodynamic potential u is regarded as a function of the two parameters {Mt, q). 
The terms 6ujd and 6u!p are defined in Eq. (p^ and (!25|) . and the effects of current 
mass m is described by the second brackets of Eq. (l42!) . Because the ground state 
is described by the two parameters {Mt,q) which minimize u, the term rn^/AG in 
the second brackets of Eq. (H2]) can be regarded as a constant. We then need to see 
remaining two terms, —{m/2G)Mt and {6u}-d + Sup). First, the term —{m/2G)Mt in 
the second brackets of Eq. fj42|) depends on Mt, not on q. The larger Mt is, the more 
markedly this term —[m/2G)Mt lowers u independently of q. Therefore, the term 
— {m/2G)Mt has possibility to make uj have minimum value at the point g 7^ 0. 

Second, let us inspect the term (^ci^d + 5uj-p)- We want to judge whether and 
5ujy tend to make uj minimized at the point g 7^ or not. To this end, the dependence 
of 5u}-c {Sui-p) on g is found by expanding Suj-q [Sujy) about g = 0, 



. 1 d\5uj^) 



2! 9g2 



■q^ + 0{q^). (43) 

q=0 
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Since the coefficient of g^, 



1 d\5u;^) 



2! 9g2 



6 V 2 / M7r2 ii 



Mt 



<7=0 



> 0, 



/A2 



(44) 



is positive, the value of duj^ at g 7^ is larger than that at g = when q is small. 
Hence, 5uj-£, tends to make uj not minimized at the point g 7^ if g is small. Similarly, 
^ti^F is expanded about g = 0, 



1 d^{6uo^) 



2! ag2 



q' + 0{q'). 



(45) 



9=0 



Since the coefficient of g^ 
1 d''{5ujY) 



2! 9g2 







9=0 





S'p 1 



(2vr) 



<0, 



(46) 



is negative, the value of 5uj-p at g 7^ is smaller than that at g = when g is small. 
Hence, 5uj-p tends to make uj minimized at the point g 7^ if g is small. Now, how 
about the sum, {5uj-q + 5ujy) ? Expanding [Sojjy + Suj-p) about g = 0, we obtain 



where the coefficient of g^ is 
1(1) 2mN,NfMt 



1 d'^{5ujD + 5ujy) 
2! 



■q^ + 0{q'), 



(47) 



q=0 



dr^e 



J1/A2 T 



d^p 1 
(27r)3 7 



0(/i 



(4J 



The equation enclosed in the above brackets is the same with the equation in the 
brackets in Eq. flB.ip . which equation takes positive value when /i < /ic(~ 0.7A), as is 
discussed in Appendix B. The coefficient Eq. (H8|) is positive if /i < /Uc(~ 0.7A). We can 
say that when g is small the value of {Sud + 6ujf) at g 7^ is larger than that at g = 0. 
Furthermore, we conffim by numerical calculations that {6ujy) + 6ujp) takes positive 
value in the range of < g < 2/i when /i < 0.7A. Eventually, we find that {6ujj) + 6ujp) 
tends to make u not minimized at the point g 7^ 0. In other words, (5c(Jd + 6ujp) tends 
not to realize the inhomogeneous chiral condensate (g 7^ 0). 

Our argument in this subsection is summarized as follows. The influence of the 
current mass m on the thermodynamic potential u is represented as —{m/2G)Mt + 
{6ujy) + 6ujf). The term —{m/2G)Mt tends to make u minimized at the point g 7^ 0. 
On the other hand, the term (^cud + ^^^f) tends to make u not minimized at the point 
g 7^ 0. Whether the introduction of the finite current mass tends to materialize the 
inhomogeneous chiral condensate (g 7^ 0) or not depends on the competition between 
the term —{m/2G)Mt and the term {Sud + Sui-p)- 
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4 Spatially inhomogeneous baryon density 



In the previous subsection, we consider the effect of the finite current mass on the 
ground state with the spatially inhomogeneous chiral condensate. The influence of 
the current mass on the thermodynamic potential is represented as —{m/2G)Mt + 
{5uj-D + 5uj-p)- There is, however, another influence which the current mass has on the 
system. In the chiral limit, the baryon density is spatially homogeneous even if the 
spatially inhomogeneous chiral condensate occurs [HlEj. In contrast, when the current 
mass is finite, the baryon density becomes spatially inhomogeneous if the spatially 
inhomogeneous chiral condensate occurs. This will be shown as follows. 

The quark number (one-thirds of the baryon number) density p can be written by 
the field ip' defined in section 2, 



P 



—i tr 



7°S'(a;, x) 




Oo ( \ + 757-3) c ^ , ^ 
T Sq{j)) Soij) + q) 



7°^o(p)^^^5o(p-g) 



(49) 



where S is the propagator of ip' , and the perturbative approximation Eq. (|T3|) has been 
used. Since p depends on r, the baryon number density is spatially inhomogeneous. The 
wave number representing inhomogeneousness of baryon number density coincides with 
that of chiral condensate up to the order 0(Kn). In the chiral limit, the baryon number 
density becomes spatially homogeneous, indeed. Why does the baryon number density 
become spatially inhomogeneous when the current mass is finite? The path integral 
representation Eq.(l7]) will help us to see this. The propagator of the field ip depends 
on the space coordinates r explicitly. In the chiral limit, however, the propagator S of 
the field ip' introduced in Eq.(I7j) does not depend on r. One can write p in terms of the 
propagator S in the mean-field approximation. On the other hand, when the current 
mass is finite, the propagator S does depend on r explicitly. 
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5 Conclusion 



We have studied how the introduction of finite current quark mass affects the ground 
state with the spatially inhomogeneous chiral condensate, Eq.(jl]), which is actually re- 
alized in the NJL model at finite density in the chiral limit. Our numerical calculations 
show that, even if the finite current quark mass is introduced, the spatially inhomoge- 
neous chiral condensate can take place. If GA^ > 6.94 is satisfied, the spatially inhomo- 
geneous chiral condensate occurs (i.e., g 7^ 0) at high density. When 5.6 < GA^ < 6.94, 
the spatially inhomogeneous chiral condensate does not occur (i.e., q = 0). The ground 
state was determined by obtaining the values of two parameters {Mt, q) which minimize 
the thermodynamic potential uo. When the current mass m takes a finite small value, 
the correction due to m was added to uj. That correction was divided into the following 
two parts. The first part is the term, —{m/2G)Mt, which depends on the constituent 
mass Mt but not on the wave number q. This term lowers u more markedly for the 
larger Mt, independently of q. Hence, this first part, — (m/2G)Mf, has possibility to 
tend to realize the spatially inhomogeneous chiral condensate. The second part is the 
term written as [Sujd + 6u}p) , which depends on both Mt and q. According to numerical 
calculations, the value of {Sujb + Sujp) at g 7^ is larger than that at g = for each 
value of Mt. Therefore, the second part, {6ljy) + 6ujf) tends not to realize the spatially 
inhomogeneous chiral condensate. Whether the introduction of the finite current mass 
tends to realize the spatially inhomogeneous chiral condensate or not depends on the 
competition between the term —{m/2G)Mt and the term (^cud + Sup). 

The introduction of the current mass infiuences not only the possibility of the spa- 
tially inhomogeneous chiral condensate but also space dependence of the baryon number 
density. In the chiral limit, the baryon number density is spatially homogeneous even 
if the spatially inhomogeneous chiral condensate occurs O E]- In contrast, when the 
current mass takes a finite value, the baryon number density becomes spatially inho- 
mogeneous if the spatially inhomogeneous chiral condensate is realized. In addition, 
the wave number vector of the spatially inhomogeneous baryon number density coin- 
cides with that of the spatially inhomogeneous chiral condensate when the value of the 
current mass is small. 

Here we give consideration to the ansatz Eq. (jl]) that describes the spatially inhomo- 
geneous dual chiral condensate. In the chiral limit, the Lagrangian is chiral invariant 
and this ansatz seems to be appropriate because (ipip) and {ipi 75 rs ip) are put on a 
chiral circle, (ipip)'^ + ('^/'^ 75 T3 t/^)^ = C^. On the other hand, when the current mass is 
finite, the Lagrangian is not chiral invariant and there might exists more appropriate 
ansatz than Eq.(jl]). In such ansatz, (ipip) and {ipi'-^^r^ip) are not put on the chiral 
circle and would satisfy more complicated relation. In the present paper dealing with 
the finite current mass, however, we have assumed the ansatz Eq.(jl]) for simplicity. 
Within the limits of this ansatz, we found the ground state with Mt > and g 7^ 0. 
Also, the baryon number density was found to be spatially inhomogeneous as well as 
the chiral condensate. If we can find more appropriate ansatz than Eq.(jl]), we would be 
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able to have thermodynamic potential whose minimum value is smaller than the mini- 
mum value obtained by use of the ansatz Eq.(jll). With such more appropriate ansatz 
than Eq.(jl]), physical quantities such as the baryon number density will be affected and 
the way of the phase transition might be influenced. Besides our ansatz, other type of 
ansatz was studied in Ref [HI [20] as discussed in section 1. In order to determine which 
ansatz is to be realized, one needs to compare the thermodynamic potentials between 
them. 

Several problems are still remained. First, it will be interesting to carry out nu- 
merical calculations with finite temperature T > 0. The phase diagram in the T — n 
plane should be affected when the current quark mass is taken to be non-zero [T5[ [20] . 
Second, the study with the inclusion of color superconductivity. In the chiral limit, 
the phase diagram in the T — fi plane including a second color superconducting phase 
(2SC), uniform chiral and non-uniform chiral phase is studied [21j. The introduction of 
non-zero current quark mass should influence the diagram. Finally, in the NJL model 
with quark mass term there appears the Nambu-Goldstone boson (NG boson) when 
the ground state with spatially inhomogeneous chiral condensate is realized. Since this 
inhomogeneous chiral condensate breaks space translational symmetry spontaneously, 
the NG boson appears [22j- In the low energy phenomena, the NG boson would play 
an important role. 



Acknowledgments 

The author thanks the Yukawa Institute for Theoretical Physics at Kyoto University. 
Discussions during the YITP workshop YITP-W-08-09 on "Thermal Quantum Field 
Theories and Their Applications" were useful to complete this work. 



16 



Appendix 



A Input parameters 

We explain in this appendix a way to set the values of the input parameters (G, m, A) 
of the NJL model in the proper-time regularization (PTR) scheme. In the vacuum 
state n = 0, the pion decay constant in the PTR scheme is represented as jl] 

and the pion mass m^r satisfies the following relation 

mlf^=m^M,. (A.2) 
The gap equation in the PTR scheme becomes [3] 

Mt = m + ^GN.NjMt F^iM^, -^), (A.3) 
where the function F2 is defined as 

1 /"oo 1 „ /-oo 1 



F2(M/, -Tx) = / f/s - e~^^^* ^ = AM dt- e'\ (A.4) 

Now, we give an arbitrary value to A firstly. Then, m and G are determined so as to 
reproduce the values /,r = 92.4 MeV and m^r = 135 MeV. When the value of A is given, 
one can calculate the value of Mt through Eq. flA.ip . Regarding the equations Eq. flA.2p 



and Eq. flA.Sp as simultaneous equations for G and m, one can obtain both the values 



of G and m. Thus, we can determine the values of G and m for a given value of A. 



B Chiral symmetry restoration and 
proper-time regularization 



In this appendix, we look into the reason why the Mt becomes less than the current 
mass, Mt < m for g = and large chemical potential /i ~ 0.7A in the numerical 
calculations. At the minimum point of the thermodynamic potential uj with g = and 
T = 0, the condition, 5uj/5Mt = 0, is held. 



Mt 



m 



AGN^NfMt 



02/ '^"5" 
Stt"' 7i/a2 s 



■ exp( 



-Mfs) 



d?p 1 



(2vr)^ 



fB.r 



where e = yp^ + Mf. This is nothing but the gap equation and it should be noted 
that the PTR scheme has been used in this equation. The first (second) term in the 
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brackets of Eq. llB.ll) represents the contribution from the Dirac sea (Fermi sea). The 
second term increases if the chemical potential n becomes larger. When the chemical 
potential reaches a certain value fic, the second term becomes equal to the first term, 

r ds-^exp{-Mfs)= [ ^,-0{f^c-e), (B.2) 

and then one has Mt = m from Eq. (lB.ip . If the chemical potential /i is larger than /i^, 
the second term becomes larger than the first term and one has Mt < m. In short the 
relationship between Mf and m in terms of size is summarized as follows. For small 
values of the chemical potential fi, Mf is larger than m. When fi reaches fic, we have 
Mt = m. If yU becomes larger than fi^, we have Mt < m. Now, let us estimate the value 
of fic assuming Mt/A <^ 1. The expansion of the left side of Eq. flB.2p about M^/A = 
takes the form |4j 

— - / ds—expi-Mfs) = — 7;\^] / dt—exp(-t) 



A^ 

87r2 



A2 

The expansion of the right side of Eq. (1B.2|) about Mt /A = is 

A2 



(B.3) 



_ 1 (Ml 

A2 j 2 I A2 



+ ■ 



47]-2 

Neglecting the order 0{Mf/K^) in Eq.(lR2]), we have 



(B.4) 



^ A ~ 0.7A. (B.5) 

Thus, we can understand roughly why Mt becomes less than m for ~ 0.7A in the 
PTR scheme. 

It would be helpful to examine the relationship between Mt and m in a different 
regularization scheme, say the three-momentum cut-off scheme In this scheme, the 
first term in the brackets of the gap equation Eq. (lB.ip is replaced by 

d?p 1 

■ d ( A3_momcntum - |p| )• (B.6) 



(27r)3 e 

The value of the chemical potential at which the equation enclosed in brackets of the 
gap equation becomes zero is 

(/^c) 3— momentum ~ -^3— momentum- (B-7) 

Therefore, Mt does not become less than m for yU ~ As-momentum in the three-momentum 
cut-off scheme. 
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